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Abstract. We discuss different aspects of the quark mass dependence of the 
chiral and deconfincment transitions. We make predictions for locating the 
deconfining critical point, by analyzing the mq dependence of the ratio of 
imaginary and real mass of the Polyakov loop. We then present an effective 
theory that allows the extraction of coupling constants from available lattice data 
for the temperature dependence of the chiral and Polyakov loop condensates at 
different quark masses. 

Theoretical developments accompanied by recent numerical results from the 
lattice have advanced our understanding of QCD. There are still numerous important 
questions however that need to be answered. For instance: Where does the QCD point 
lie in the (ms,mg) and in the (T,mq) phase diagrams? Which mechanism drives the 
transition for a given quark mass? Is the QCD transition more deconfining or chiral 
symmetry restoring? See ^ for related discussion. As a step towards answering these 
questions, first, we discuss the location of the deconfining critical point in the (T,mq) 
plane. Then, we introduce a general theory for the Polyakov loop and chiral field at 
finite quark mass, and discuss its fit to lattice data. 

The pure gauge sector of QCD with N colors, i.e. niq = oo, has global Zjv 
symmetry associated with the center of the SU{N) gauge group. The Polyakov loop 
£, charged under Z^, serves as order parameter for the deconfining transition. £ is 
real for N = 2 and the Z2 symmetry breaking deconfinement transition is of second 
order. For N — 3 the Polyakov loop is complex and transforms as £ — > exp(27ri/3)€. 
Accordingly, the expectation value {£) = at low temperatures, when Z3 is unbroken, 
and {£) 7^ above the deconfinement critical temperature Td . The deconfinement 
transition for iV = 3 is weakly first order. This has been verified by lattice QCD jj. 

In an effective theory for the Polyakov loops, the Polyakov Loop Model (PLM) 
has been constructed!. The potential for iV = 3, up to quartic terms, is 

no.-.,^-|i.p-|^i±^:! + l(m^ (1) 

In a mean field analysis all coupling constants, except 62 (T'), are temperature 
independent. There are two masses, associated with the real £r = Re £ and imaginary 
£i = Im £ parts of the Polyakov loop. These are defined from the two-point functions 

(4(x)4(0)) -£l- ^M-mrx) ^ (^^(^)^^(o)) ^ ^MzH^ (2) 
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and thus characterize the exponential tail of the correlators. {£) = £o > is real. 
Denote the two minima i" = ^(ly) and ^+ = 1{T^). The masses at are then 

ml = -62 - 263^+ + 3{e+y , m\ = -62 - 263^+ + (^+)' . (3) 

For &i = the potential is Z3 symmetric. The solutions are 

^0=0, ^ = ^^3, h^ = -\hl. (4) 

and the ratio of the masses at T'^ is nii/nir = 30]. This value is twice the one 
expected from lowest order perturbative calculations, valid at high temperatures. The 
PLM thus predicts an increase of mi/rrir as T ^ T^, prediction supported by the 
lattice In principle, mi/rrir can receive contributions from pentic and hexatic 
terms (see discussion in 0]). 

What happens to rrii/mr when dynamical quarks are added in the theory? How 
does this ratio change for finite quark masses? For any finite the symmetry 
is explicitly broken, and the linear term bi is turned on in This corresponds to 
tilting the potential and shifting the degenerate minima at Td towards each other, and 
thus the transition becomes weaker first order. An increasing symmetry breaking is 
equivalent to decreasing the quark masses from infinity. An explicit functional relation 
for bi{mq) fitted for lattice data, together with the mq dependence of Td is discussed 
in For some value of mq there is only one minima and the transition is of second 
order. This is the deconfining critical point D, shown in the left panel of Fig. ^ 
When further decreasing the quark mass the phase transition becomes a crossover. To 
determine the infiuence of quark mass we derived a set of differential equations 
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dhi m^V^95iV' a6i miy^dbi^)' dbi 

and solve them together with the initial conditions Q numerically. Here V"{£~) = 
and V"{£~^) = m^, and 63 = 0.9. Note that the third term is always negative, leading 
to decreasing Td with mq (see the slope of the upper 1st order critical line in the left 
panel of Fig. The mass ratio at the critical temperature in terms of the symmetry 
breaking is shown in the right panel of Fig.^ We find that as the quark mass decreases 
from infinity the ratio slightly increases. The second order deconfining critical point D 
is reached for 61 = 0.027 and here the mass ratio becomes divergent. This means that 
only the real part ir , not the ii , becomes massless at D, and thus Z3 is a symmetry 
group here. Following the behavior of mi/mr on the lattice would allow for locating 
D, and also to determine the quark mass that this point corresponds to. 

Another well studied sector of QCD is in the limit of zero and small quark masses. 
In this regime Zn is always broken, but chiral symmetry is restored above a critical 
temperature Tc- The order parameter is the chiral condensate a. The transition is 
second order for two, and first order for three fiavors of massless quarks. For finite 
quark masses chiral symmetry is explicitly broken. Increasing mq from zero means 
weakening the first order phase transition. For some mq the first order critical line 
turns into a second order chiral critical point, C in the left panel of Fig. For even 
larger mq the transition is a crossover. Effective field field theories have been used to 
analyze this sector 0. 

Deconfinement and chiral symmetry restoration are two phenomena very different 
in nature. For realistic quark masses both the chiral and Zjy symmetries are explicitly 
broken, and likely that both transitions are in the crossover domain [7]. Here neither 
fT, nor £ can play the role of a true order parameter. However, lattice simulations 
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show that the susceptibihties associated with these two quantities peak at the same 
temperature, when quarks are in the fundamental representation of the gauge group 
1^, indicating that the transitions occur at the same critical temperature. Similar 
results were obtained also in terms of the chemical potential TTF . 

In we provided an explanation within a generalized Ginzburg-Landau theory 
to why for rriq = chiral symmetry restoration leads to deconfinement. The analysis 
was based on the following general idea: The behavior of an order parameter induces a 
change in the behavior of non-order parameters at the transition |12j , via the presence 
of a possible couphng between the fields, giia^. In JT] we assumed gi > 0, which is 
confirmed by extracting it from the lattice |13|. 

Here we extend the analysis of JJj for finite quark masses. We can then study not 
only the qualitative relation between deconfinement and chiral symmetry restoration 
in terms of quark masses, but also can provide quantitative statements after extracting 
the coupling constants from lattice data 13^. Recent lattice results for the temperature 
dependence of the chiral condensate for three degenerate massive quarks were reported 
in 17 . The first results for the behavior of the renormalized Polyakov loop as 
a function of temperature, under the same lattice conditions, were reported in 

We now write down the most general renormalizable effective Lagrangian that 
can be built from the chiral field cr and the Polyakov loop £. The contribution to the 
potential V — Vpi + Vch + Vint from the Polyakov loop is 

VMl)=9oi-^ji'-^-^i' + le'. (6) 

Here go = —bi §. For the chiral potential we adopt the linear a model, discarding all 
fields but the a, since in relation to our discussion these will not play a role |13) : 

= + (7) 

Here H — f-nrn^ , with TOtt and the vacuum pion mass and decay constant. The 
possible interactions terms are 

§ This potential is written only for the real part of the Polyakov loop, since we are interested in the 
behavior of its expectation value, which can always be chosen real at /i = 0. 
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In a mean field analysis the field expectation values to first order in rrig are 

{e)c.-J^-^{ar-^H{a) (9) 
and (a) is the solution of 

For if = we recover the results of pl|, as expected. For light quarks of three 
degenerate flavors, the Polyakov loop mass is identified as the inverse of the screening 
radius and is extracted from the lattice ^H]. For the a mass we assume m^{T) = 
a {bT/Tc — 1) as is customary in a Ginzburg-Landau theory. Then using and (|10|l 
we fit the data from (16j and ^1]. The couplings extracted this way confirm the 
assumption made in in that gi > and go < 0. Accordingly, for massless quarks 
(H = 0) the decrease of the order parameter, the chiral condensate induces an increase 
in the Polyakov loop, as one can see directly in Which field drives the transition in 
the case of massive quarks can be infered from the knowledge of the couplings. These 
results will be presented elsewhere |13j . 

With our effective model, we will be able to make predictions relevant for the 
phenomenology of heavy ion collisions. The extension of this analysis to the chemical 
potential axis of the phase diagram is relevant not only in the context of the critical 
point in the (T, ^) plane, but also in the understanding of the phase structure of QCD. 
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